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Abstract 

In this paper (the second of a series) we extend our calculation of a 
classical fixed point action for lattice SU(3) pure gauge theory to include 
gauge configurations with large fluctuations. The action is parameterized 
in terms of closed loops of link variables. We construct a few-parameter 
approximation to the classical FP action which is valid for short correla- 
tion lengths. We perfor m a scaling test of the action by computing the 
quantity G = LyJ <r(L) where the string tension er(L) is measured from 
the torelon mass /i = La(L). We measure G on lattices of fixed physical 
volume and varying lattice spacing a (which we define through the decon- 
finement temperature). While the Wilson action shows scaling violations 
of about ten per cent, the approximate fixed point action scales within 
the statistical errors for 1/2 > aT c > 1/6. Similar behaviour is found for 
the potential measured in a fixed physical volume. 
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1 Introduction and summary 



This paper is the second in a series M describing the construction and testing of 
a fixed point (FP) action for lattice ST/ (3) pure gauge theory. In I we derived the 
general FP equations for the action and for selected operators. We discussed 
their properties on the classical level and in 1-loop perturbation theory. We 
solved the FP equation for the action and the Polyakov loop on fields with 
small fluctuations. The FP action is a perfect classical action on the lattice 
and, as we argued in I, it is perfect even in 1-loop perturbation theory. It 
is the [3 — oo point of the renormalized trajectory and, as such, will be the 
basic ingredient of future attempts to construct perfect actions. This should 
be the ultimate goal in order to reach the continuum limit without losing the 
battle against critical slowing down, memory problems and other theoretical or 
computational issues. 

The aim of this paper is more modest. It is natural to expect, and this expec- 
tation is supported by a numerical study of the non-linear cr-model Q , that the 
FP action works much better than the Wilson action even at small correlation 
lengths. In this paper we perform a scaling test and make a comparison. 

Scaling means that all physical dimensional quantities show the same func- 
tional dependence on the gauge coupling. It is different from asymptotic scaling 
where in addition to scaling we require that this functional dependence be de- 
scribed by the 2-loop A parameter. Asymptotic scaling properties can be signifi- 
cantly improved by a suitable non-perturbative redefinition of the bare coupling 
[[|, H . In this paper we are not concerned with asymptotic scaling. Scaling prop- 
erties can only be changed by modifying the lattice action. A redefinition of the 
coupling constant used in the lattice action during simulations will also shift 
parameters in the lattice action away from their naive free-field values. Sim- 
ulations done with these new actions may show better scaling properties. An 
example of the use of such actions is the self-consistent tadpole-improvement 
program for nonrelativistic QCD which has been applied with great success to 
studies of charm and bottom quark systems Q. 

Scaling violations can be different for different quantities. The ratio of cer- 
tain observables might show scaling at smaller correlation length than others. 
Scaling requires the universal behavior of all observables. Of course, we cannot 
test all observables, but we can ask whether observables which show scaling 
violations when computed using the Wilson action scale better when using a 
FP action. 

The first step is to derive a parametrization for the FP action which is 
sufficiently simple and approximates well the FP action on rough configurations 
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which typically enter our simulations. This is a multistep procedure which we 
discuss in some detail in Sections 2 and 3. At the end we obtained an 8 and a 16 
parameter approximation of the FP actions for RG transformation of type I and 
type II, respectively. (These RG transformations were defined and discussed in 
detail in I.) These parametrizations are constructed in terms of two loops, the 
plaquette and the twisted perimeter-six loop. Powers of the real and the square 
of the imaginary part of the loop-trace enter the action with coefficients c\, c 2 , ... 
and d\, d 2 , respectively as given in Tables | and |. 



Table 1: Couplings of the few-parameter FP action for the RG transformation 
of type I. The couplings d are zero. 



operator 


Cl 




C2 




C3 c 4 


Cplaq 
^6—link 


.523 
.0597 




.0021 
.0054 




.0053 .0167 
.0051 -.0006 


Table 2: Couplings of the few-parameter FP action for the RG transformation 
of type II. 


operator 


Cl 


C2 




C3 


C4 c 5 


Cplaq 
^■6 — link 


.6092 
.0489 


.0478 
-.0146 




-.0470 
.0393 


.0295 -0.0038 
-.0216 .0039 


operator 


di 


d 2 




d 3 




Cplaq 
^6 — link 


.0921 
.0042 


-.1487 
-.0034 




.0397 
.0042 





We pushed the scaling test to the extreme, going to configurations with very 
large fluctuations (corresponding to Wilson (3 = 5.1). This extreme situation 
might be interesting for finding the connection with strong coupling expansions 
H but is, presumably, not practical in standard applications. On very rough 
configurations it is increasingly difficult to find a simple parametrization and, 
in addition, the correlation length is so small that the signal disappears rapidly. 

In the scaling test which we performed with the 8 parameter action in Table 
[l], the critical temperature T c was used to set the physical scale. We determined 
the critical coupling constant /3 c (N t ) for for N t = 2, 3, 4 and 6 (Table||) and fixed 
the lattice spacing at these coupling values as a = l/(T c N t ). The corresponding 
critical couplings for the Wilson action are known to good precision [Q . 

Next we measured physical observables at couplings (3 c (N t ) in fixed, finite 
physical volumes. This way we can avoid infinite volume extrapolations. For 
one scaling test we consider the quantity G = L-*/ ct(L) where a(L) is the string 
tension on an L 3 volume computed from the exponential fall-off of the Polyakov 
line correlator (torelon mass). Taking L = rN t where r is some conveniently 
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Figure 1: Scaling test for the Wilson action (crosses) and FP action (squares); T c 
is defined in infinite spatial volumes. The diamonds are from zero-temperature 
simulations by Ref. 8 and the diamond and the star at a = are extrapolations. 



chosen aspect ratio, the physical volume is kept fixed at V = (r/T c ) 3 . The ratio 
G is independent of the bare coupling (resolution) in the scaling (continuum) 
limit. Any variation of G is due to lattice artifacts. Figure |l| shows G(L) with 
r — 2 as a function of aT c for the Wilson and for the 8 parameter FP action. The 
inner and outer error bars on the FP points show their statistical uncertainty, 
and the combined uncertainty from statistics and in f3 c . The Wilson action 
shows a scaling violation of about 10% between N t — 2 and 6. No scaling 
violation above the statistical errors is seen for the FP action. 

No simulations using the Wilson action at small lattice spacing precisely 
match our scaling test using the torelon mass on fixed physical volumes. The 
ones which come closest are the very recent ones of Boyd. et. al. H. These 
authors present measurements of the quantity ^/a/T c at the critical couplings 
for Nt = 4, 6, 8, and 12, for which the string tension has been computed via 
fits to Wilson loops on large (zero temperature) lattices. (They also present 
new measurements of critical couplings at Nt — 8 and 12.) The string tension 
from Wilson loops is an upper bound on cr{L); a(L) is reduced from a by the 
zero-point fluctuation term which for large L in a string model is || 

<L) = a-^ + .... (1) 
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The diamonds in Fig. [I] show 2yfa/T c from Ref. ||. The authors of that 
work present an extrapolation to a = 0, which we also display along with the 
a = limit of Ly^a(L) at r = 2 calculated from Eq. ([!]). 

Torelon masses at small lattice spacing have been presented in the literature, 
but they are not at couplings which are critical couplings for deconfinement for 
integer N t 's. Our attempts to interpolate them into scaling tests which match 
our simulations generally resulted in uselessly large uncertainties, and we do not 
show them here. 

Results are similar from Fig. || where the aspect ratio r = 3/2 data are 
plotted. In both cases the question remains unanswered, whether the constant 
observed with the FP action at very low resolutions really agrees with the con- 
tinuum value. 
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Figure 2: Scaling test for the Wilson action (crosses) and FP action (squares) on 
aspect ratio 3/2 lattices. T c is defined in infinite spatial volumes. The diamond 
is a Wilson action simulation by us at (3 — 6.0 and the star is an extrapolation 
using two-loop scaling to the correct (Ref. 8) j3 c — 6.06. 

Another scaling test is offered by the q — q potential V(r; L) in a volume L 3 
where L is chosen here to be L = 2N t . Unfortunately, this quantity is contam- 
inated by cut-off effects for r ~ a coming from non-perfect sources. The FP 
Polyakov loop has not been constructed yet beyond the quadratic approxima- 
tion (paper I), so we measured simple Polyakov loop correlations. The different 
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behaviour of the potentials obtained with the Wilson and FP actions is well 
demonstrated by Figs. || and ||, nevertheless. 

4 | — , — , — , — , — | — , — , — , — , — | — , — , — , — , — | — , — , — , — , — | 
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Figure 3: Potential V(r)/T vs. rT for the Wilson action at (3 C {N T = 2) (filled 
circles) and (3 c (Nt = 4) (crosses). 

In the potential of the FP action the rotational symmetry violation at large r 
drops below the statistical errors even at j3 c (N t ) = 2 but remains clearly present 
for the Wilson action. (Compare, for example, the predictions at r = 1.5/T C 
coming from r = (0,0,3) and f — (1,2,2) at /3 c (Nt) = 2.) The non-physical 
additive constant in the potential was fixed by the convention V(r = l/T c ) = 0. 
The lack of scaling for the Wilson action is evident, while the long distance 
behaviour of the potential is consistent with scaling for the 8 parameter FP 
action. 

The outline of the paper is as follows: In Section 2 we describe the con- 
struction of a classical FP action which is appropriate for gauge configurations 
whose fluctuations are not small. In Section 3 we give an approximate simple 
parametrization of the FP action which represents the FP on gauge field config- 
urations which arc typical at the correlation lengths of our simulation. Finally, 
in Section 4 we give some details of the scaling tests using this action. 
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Figure 4: Potential V{r)/T vs. rT for the FP action at (3 C {N T = 2) = 2.975 
(filled circles) and (3 c {Nt = 4) = 3.52 (crosses). 



2 Numerical analysis of the FP equation and the 
problem of parametrization 

We consider an SU(N) pure gauge theory^ on the lattice and consider the RG 
transformation 



e 



WW = / DU exp {-13 (S(U) + T(U, V))} , (2) 



where U is the original link variable, V is the blocked link variable and T(U, V) 
is the blocking kernel that defines the transformation. At (3 = oo the transfor- 
mation becomes a steepest decent relation 

S FP (V) = mm (S FP (U) + T(U,V)) , (3) 

where 

T(U, V) = -^Y1 [ ReTr (V»MQl(n B )) - max (ReTr(^gt)}] . (4) 



3 The actual numerical analysis and simulations in this paper were done for SU(3). The 
equations are written for general N if not indicated otherwise. 
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In Eq. (g) W e SU(iV) and the N x N complex matrix Q^ub) is the block 
average. Its form depends on the block transformation. For the two different 
block transformations (type I and type II) which arc considered here Q is defined 
in Figs. 5 and 6 in I. For further details we refer also to I. 

In I we solved Eq. (0) up to cubic (quadratic) order in the vector potentials 
for type I (type II) block transformation. If the configuration V has large 
fluctuations one has to use numerical methods to solve the FP equation. 

For numerical simulations we need a parametrization for the FP action which 
is sufficiently simple to make the computational overhead acceptable. It turns 
out, however, that the process of solving Eq. ^ requires, as an intermediate 
step, a more general parametrization. This is the problem we discuss first. 



2.1 Parametrization of actions 

Our parametrization is based on powers of the traces of the loop products Uc = 
^cU^n) where C is an arbitrary closed path 

S(U) = jj^2( Cl (C)(N-ReTr(U c )) + 02(C) (JV - RcTr(C/ c )) 2 + ... 
c 

+d 2 (C)(ImTr([/ c )) 2 + d 4 (C)(ImTr([/ c )) 4 ...). (5) 

We also checked other possibilities. The success of a parametrization of the 
er-model and the form of exact classical solutions in both models suggest building 
gauge invariant quantities from Uc by using the angles after diagonalization, 
U d c iag = {e ie \e i9 \e l(>3 ) for N = 3 Jig. For SU(3) gauge theory this type of 
parametrization did not work better than the one in Eq. (|5|) and was, therefore, 
abandoned. We also tried to include the product of traces of different loops. 
It is quite difficult to keep track of the large number of topological possibilities 
and, in those cases we checked, the parametrization did not become better. The 
following results are based on the parametrization in Eq. (^) . 

Using optimized parameters in the block transformation the quadratic and 
cubic part of the FP action is short ranged (paper I). So, we decided to restrict 
the set of loops C to the 28 topologically different loops that are length 8 or less 
and fit into a 2 4 hypercube. They are tabulated in Table ||. Only 7 (15) of the 
28 are independent on the quadratic (cubic) level. For the RG transformation 
of type I the FP action is known up to cubic order analytically which can be 
used to fix 15 of the ci(C) coefficients by a \ 2 fit- We decided to keep only 12 
operators since the fit did not improve further. For the RG transformation of 
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type II 7 ci(C) couplings were fixed this way since the analytic cubic result is 
not available. Table [I| and Table || contain the c\{C) coefficients obtained. We 
shall denote the action where these couplings are kept only by S pp . 

Table 3: Loops used in the construction of the FP action. One particular ori- 
entation is shown; all possible orientations are included in the actual operators. 
label path 

1 x,y,-x,-y 

2 x,y,y,-x,-y,-y 



3 


x, 


yz 


-y 


-x, 


-z 








4 


x, 


yz 


,-x 


-y, 


-z 








5 


x, 


y> 


x, 


y> 


-x, 


-y> 


-x, 


-y 


G 


x, 


y, 


x, 


y> 


-x, 


-x, 


-y 


-y 


7 


x, 


y> 


x , 


z, 


-y> 


-x, 


-x, 


-z 


8 


x, 


y> 


x, 


z, 


-x, 


-z, 


-x, 


-y 


9 


x, 


y, 


x, 


z, 


-x, 


-y> 


-x, 


-z 


10 


x, 


y, 


x, 


z, 


-x, 


-x, 


-z, 


-y 


11 


x, 


y, 


x, 


-y, 


-X 


, y. 


-x, 


-y 


12 


x, 


y, 


y. 


x, 


-y, 


-x, 


-x, 


-y 


13 


x, 


y, 


y. 


z, 


■y, 


-z, 


-x, 


-y 


14 


x, 


y, 


y. 


z, 


■y. 


-y 


-x, 


-z 


15 


x, 


y, 


y. 


z, 


■y, 


-x, 


-z, 


-y 


1G 


x, 


y, 


y. 


z, 


•x, 


-y, 


-y 


-z 


17 


x, 


y, 


y. 


-x, 


z, 


-y, 


-y 


-z 


IS 


x, 


y, 


y. 


-x, 


-X 


, -y 


-y 


x 


19 


x, 


y, 


z, 


t, - 


z , 


-t, 


-x, 


-y 


20 


x, 


y, 


z, 


t, - 


z, 


-y, 


-x, 


-t 


21 


x, 


y, 


z, 


t, - 


z, 


-x, 


-t, 


-y 


22 


x, 


y, 


z, 


t, - 


z, 


-x, 


-y 


-t 


23 


x, 


y, 


z, 


t, - 


y> 


-x, 


-t, 


-z 


24 


x, 


y, 


z, 


t, ■ 


y> 


-x, 


-z, 


-t 


25 


x, 


y, 


z, 


t, ■ 


x, 


-y, 


-z, 


-t 


2G 


x, 


y, 


z, 


-y, 


-z 


y, 


-x, 


-y 


27 


x, 


y, 


z, 


-y> 


-X 


y, 


-z, 


-y 


28 


x, 


y, 


-X 


, -> 


, X 


, y, 


-x, 


-y 



2.2 Numerical minimization 

To go beyond leading order we solved the fixed point equation Eq. (p|) nu- 
merically. The procedure is the following: First we create an arbitrary SU(3) 
configuration {V}. That serves as a coarse configuration. Our goal is to find a 
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Table 4 


: Couplings 


of the fixed point action of Type I given in 


terms of loops 


defined 


in Table 3. 


The quadratic coupling 


s are labelled by C\ 


, and the coef- 


ficients of higher powers of (3 — TW(C) fro 


m the fit described 


in Sec. 2.2 are 


also shown. 








loop 


Cl 


C2 


C3 


c 4 


1 


.6744 


.1067 


-.1072 


.0186 


2 


-.02 


.0424 


-.0304 


.0055 


3 


.012 


.0361 


-.0228 


.00436 


4 


.005 


.1132 


-.0632 


.00937 


5 


-.0031 


.0765 


-.0562 


.0109 


8 


-.0035 


.00525 


-.00418 


.00076 


9 


.0027 


.00402 


-.00642 


.00130 


15 


-.0024 


-.0356 


.0187 


-.00268 


16 


.0013 


-.0023 


-.00373 


.00132 


20 


.003 


.0193 


-.0119 


.00211 


23 


.0032 


.0187 


-.0128 


.0024 


25 


.0024 


.0126 


-.0115 


.0021 



fine configuration {U} according to Eq. (Q). Since we do not know S FP (and 
even if we knew it, it presumably would be too complicated to use in a numer- 
ical minimization), we approximate S FP with a simpler action Sq. We choose 
Sq such that it is sufficiently simple for minimization but close to the fixed 
point action. The minimization of Eq. (||) gives the fine configuration {Uo} as 
a minimum 

S'(V) = min(5 (C/) + T(U, V)) = S Q (U ) + T(U , V). (6) 

If So is close to the fixed point action S FP , than Uo is close to {U FP } and 
up to quadratic corrections 

S FP (V)= S FP {Uo)+T{Uo,V) 

= S FP (Uo) + (S'(V)-So(Uo)). (7) 

This equation can be used to calculate the value of the fixed point action on the 
course configuration {V} if S FP is known on the fine configuration {Uo}. 

One might think we got nowhere this way but this is not so. One observes 
that the minimizing configuration U FP (V) (which is close to Uo(V)) has much 
smaller fluctuations than the configuration V. As Eq. (^) shows, T(U FP (V), V) 
is non-negative and so Eq. (§) implies that S FP {V) > S FP {U FP {V)). From 
this inequality it follows that the action density on the fine lattice is at least 
1/2 4 times smaller than on the coarse one. Actually, the observed ratio is even 
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Table 5: Couplings of the fixed point action of Type II, labeled as in Table 4 . 



loop 


Cl 


C2 


C3 


C4 


1 


.5454 


-.0837 


0.0459 


-.0076 


2 


.0094 


-.0335 


.0161 


-.0026 


3 


.0 


.0639 


-.0403 


.0076 


4 


.0419 


.0594 


-.0485 


.0091 


5 


.0 


.0088 


-.0060 


.0010 


8 


.0 


-.0022 


.0008 


-.0001 


9 


.0019 


.0033 


-.0019 


.0003 


17 


-.0011 


.0035 


-.0015 


.0002 


18 


-.0041 


-.0081 


.0040 


-.0006 


20 


.0 


-.0109 


.0074 


-.0011 


23 


.0 


-.0213 


.0111 


-.0014 


25 


.0043 


.0054 


-.0039 


.0007 



smaller, it lies between 1/30 and 1/50 on the configurations we considered^. On 
the configurations Uq(V) which have very small fluctuations the FP action can 
be replaced by S FP defined at the end of Section 2.1 [J. Eq. (Q) gives then the 
value of the FP action on the configurations V. 

The procedure described above could fail if So in the original minimization 
was not close to the FP action. We found that by carefully adjusting the cou- 
plings of Sq one can make the difference S fp (Uq) — So(Uo), that measures this 
correction, small. After a certain amount of trial and error we were able to 
reduce these "perturbative" corrections to less than 0.4% and 0.3% of the ac- 
tion for the type I and type II RG transformations, respectively. In fig. |^ we 
illustrate this fact for the type I RGT. 

In the case of the RG transformation of type I we generated about 400 
configurations on 2 4 and 4 4 coarse lattices using the Wilson action, with coupling 
(3 ranging from 5.0 to 7.0. For type II we created more than 700 configurations 
in the coupling constant range 5.0 to 50.0. Following the steps of the discussion 
above we associated to each configuration a real number S FP (V), i.e. the value 
of the FP action on the configuration. 

Our next task is to find a parametrization for S FP which reproduces the 
action values to an acceptable accuracy and, at the same time, is sufficiently 
simple to be used in numerical simulations. This problem will be discussed in 

4 In order to appreciate how small is this ratio let us remark that the Wilson action density 
(3 — ReTrC/p) is reduced by a factor of ~ 2 as /3 is changed from 5.0 to 10.0. 

5 One might and we did include further terms in order to cope with the small corrections 
which arise when V is very coarse and so the fluctuations of Uo(V) are not so small anymore. 
In order to keep the discussion straight we do not follow this point further. 
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Figure 5: The perturbative corrections 5Spt = S FP (Uo) — So(Uo) for a set of 
configurations using the type I RGT. 



Section 3. Beyond that we would like to make the information stored in the 
generated configurations and action values available to those who might want 
to create their own simple paramctrization for simulations. For this purpose 
we parametrized S FP as in Eq. (^) keeping a large number of loops from the 
set in Table |[ The couplings are presented in Table || and Table || for the RG 
transformations of type I and type II, respectively. Of course, these complicated 
actions are not practical to be used in simulations. They are used only to encode 
information which would be difficult to present otherwise. The actions in Table [| 
and Table || satisfy the FP equation Eq. (||) on our configurations to a reasonable 
precision. The difference between the left and right hand sides of Eq. (|J) is no 
more than 2% of the action, and in most of the cases the difference is much less. 
The quality of the fit is illustrated in fig. [| 



12 



0.4 




_Q 4 I I I I I I I I I I I I I I I I I I I I I I I I I 

4 6 8 10 12 14 

S 



Figure 6: SS (the difference between the two sides of Eq. (0)) vs. S m i n for the 
many parameter fit of Table 4 for the type I RGT. 

2.3 Technical remarks on numerical minimization 



In order to avoid programming errors two codes were written independently 
(this applies for all the essential parts of the numerical analysis presented in 
this paper). In both cases the lattice was swept through repeatedly, making 
a random local change and accepting the change if it lowered the action. The 
local change consisted a random rotation on a fine link, or on two neighbouring 
fine links at the same time. The rotations were parametrized and the action 
was interpolated quadratically in the parameters to find a set of values which 
minimized the action in this local few-parameter space. The configuration was 
overwritten only after checking that the action was really lowered. 

In order to speed up the minimization process it was useful to begin with a 
fine configuration close to the final answer. One strategy for constructing such a 
configuration was to drop the JieTxWQ^ normalizing term from T(U, V) during 
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the minimization and to relax for 20 to 50 sweeps. This typically brought the 
action to within ten per cent of its minimum value. Then 20 to 50 sweeps of 
the full minimization were required to achieve convergence at the fourth digit of 
accuracy. Another method was to build a fine starting configuration using the 
connecting matrix Z (Eq.(25) in I) valid on configurations with small fluctua- 
tions. This method, which required temporary gauge fixing, worked suprisingly 
well even on rough configurations. 



3 The few-parameter FP action 

For numerical simulation we need a parametrization which is simple and ap- 
proximates the FP action sufficiently well on those configurations which are 
typical in the actual calculation. In order to find such a parametrization we 
used coarse configurations generated by the Wilson action between j3 = 5.2 and 
6.0 and fitted S FP (V) with powers of the traces of only 2 loops. 

The two loops we used in our few-parameter fit were the plaquette and loop 
#4 in Table ^ the twisted 6 link operator (x,y,z,-x,-y,-z). For RG transformation 
of type I we kept a,i = 1, ...,4 and put di = 0. The fit for the 8 parameters 
are listed in Table [l]. For RG transformation of type II we parametrized with 
Ci,i — 1,...,5 and di,i = 1,...,3. The corresponding 16 fitted parameters are 
given in Table g. 

The accuracy of the fit is illustrated for the 8 parameter action in Fig. [7] 
where SS = S flt (V) ~ S FP (V) is plotted as a function of S FP . The quality 
of the fit deteriorates for smaller correlation lengths (larger S FP ) but remains 
within the 5% range. 

As is evident from Tables [j] and [|, the FP action is dominated by the pla- 
quette term with a small contribution of the 6 link operator in the lowest power. 
Can we neglect higher powers of loops altogether? To illustrate the importance 
of these small coefficients, in Fig | we plot SS = S flt (V) - S FP (V) calculated 
with only the quadratic c\ coefficients. It is obvious that the seemingly small 
Ci, i = 2, 3, 4 terms play a very important role in fitting the FP action. 

The parameters in Tables [l] and |^ were chosen to approximate the FP action 
on configurations typical of a correlation length range occurring in pwuson ~ 
(5.2, 6.0). If we wanted for some reason to perform simulations at 10 -3 fm, we 
would use a parameterization of the FP action which was valid at the small 
lattice spacing. 
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Figure 7: SS = S fit (V) - S FP (V) vs. S FP for the 8 parameter fit. 



4 Simulations with approximate FP actions 



The FP action is the classically perfect lattice action. As we argued in I, it is 
perfect even in 1-loop perturbation theory. These features make it plausible to 
expect that the FP action will have improved properties even at small corre- 
lation lengths. We present a detailed scaling test for the FP action of the RG 
transformation of type I using the approximate parametrization in Table [l]. 
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Figure 8: 8S vs. S m i n using only the quadratic coefficients of the action. 
4.1 Simulation algorithms 

We have adapted two algorithms for simulation. First, we have written a mixed 
Metropolis/overrelaxation algorithm which acts on SU(2) subgroups, in com- 
plete analogy with a standard Wilson action code. No special optimization was 
done. The program runs with 8 parameters about a factor of 5 slower than with 
1 parameter corresponding to the Wilson action. This latter is a factor of 1.4 
slower than a highly optimized Wilson code. There are three reasons for the loss 
of speed. The first is the extra loops needed by the new action. The second is 
that because the action is not linear in a link variable, one cannot sum the links 
in the staples ahead of time, but must do so for each new update. Finally, again 
because the action is not linear in the link variables, we have not been able to 
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invent a microcanonical overrelaxation step, and after each overrelaxation step 
we must recompute the action and perform a Metropolis accept/reject step. The 
recomputation is unnecessary for the Wilson action. We have not precomputed 
any parts of the action for greater speed; the program is designed for arbitrary 
actions. 

All our running on serial machines has used this algorithm. We typically do 
a mix of four overrelaxation updates on each of three SU{2) subgroups on each 
link, followed by a one-hit Metropolis update (this is a "sweep" as used below), 
mimicking a conventional updating pattern for the Wilson action. We have 
not done long enough runs to have measured reliable simulation autocorrelation 
times. They appear to be about five sweeps for the plaquette, measured near 
the Nt — 4 critical coupling. The autocorrelation time was negligible for the 
correlators of Polyakov loops. We typically measure all quantities once per 
sweep. Of course, the whole point of the perfect action program is to be able to 
run on small lattices with small correlation lengths where critical slowing down 
is absent. 

The overrelaxed/Metropolis algorithm is rather inconvenient for parallel ma- 
chines. The action does not parallelize on two checkerboards due to the presence 
of the six-link term, and because the action includes more than loops in the fun- 
damental representation, much temporary storage is needed. Neither problem 
is insurmountable, but we avoided them by using the Hybrid Monte Carlo algo- 
rithm pl] |. The HMC code we have written for the eight-parameter FP action 
is about a factor of seven slower than a highly optimized Wilson action HMC 
code [Oj. The difference in speed is easily understood as coming entirely from 
the many extra loop terms in the action. 



4.2 The method of testing the FP action 

When scaling holds the ratio of two physical masses is independent of the lattice 
resolution, which can be changed by tuning the coupling (3. One of the physical 
masses we took is the critical temperature T c which is used to set the physical 
scale. The lattice spacing a and T c are related in the usual way, by finding the 
critical coupling (3 c {N t ) at which the confinement-deconfmemcnt transition oc- 
curs on a lattice of size N% x N t (with N s ^> N t ). All our additional simulations 
were done at the discrete (3 values (3 c (N t ), where the lattice spacing is already 
fixed and given by a = l/(T c N t ). 

For the other physical mass we measured the torelon mass. A somewhat less 
quantitative test is offered by the potential which we also measured. 
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In the first case the scaling test is the quantity G = Ly/a(L) where a(L) 
is the string tension on an L 3 volume measured through the torelon mass. We 
performed the numerical experiment at two different aspect ratios L/Nt = 2 
and 3/2. 

In the second case we studied the potential V(r; L) in a volume L 3 , where L 
was chosen to be L = 2N t . This test is less quantitative, since, rather than using 
FP Polyakov loops which we have not yet constructed beyond the quadratic 
approximation (paper I) , the potential was extracted from simple Polyakov loop 
correlations. This is correct only for r >> a where it makes a negligible error 
that the Polyakov loops create sources whose distance is not exactly r. 



4.3 Critical couplings for deconfinement 

We found that the Columbia group's definition |l4j of (3 C is the easiest to 
use on small volumes. Specifically, we make simulations of up to several tens 
of thousands of sweeps, recording the value of the real and imaginary parts 
of the Polyakov loop averaged over the lattice. We form the angle 9 from 
tan 6 = InxP/ReP and measure the fraction of the time /20 that 9 lies within 
±20 degrees of one of the ordered orientations expected for the Polyakov loop in 
the high temperature phase. We define (3 C as the place where the deconfinement 
fraction ^ 

Id = 2/20 - g ( 8 ) 

equals fifty per cent. Most of our simulations are done on lattices of spatial 
size L between 2 or 3 times N t . On these sizes close to (3 C the lifetimes of the 
various Z 3 vacua average a few hundred sweeps. As we will see, it is necessary 
to measure (3 C to within ±0.005 to make meaningful scaling tests. We found 
alternate diagnostics such as histograms of the average Polyakov loop modulus 
just less precise than fjj measurements. 

Infinite volume (3 C values are found by extrapolating finite volume f3 c values 
linearly in 1 /volume and also by looking at crossings in a plot of Binder cumu- 
lants, which for us is B = 1 — (L 4 ) /(3(L 2 ) 2 ), where L 2 is the absolute square 
of the volume-averaged Polyakov loop on a configuration, L A = (L 2 ) 2 , and () is 
an average over a simulation run. 

We have measured critical couplings on lattices with Nt — 2, 3, 4, and 6. The 
spatial volumes and critical couplings are shown in Table ||. Examples of plots 
of deconfinement fractions are shown in Fig. |j. A plot of j3 c vs. 1/volume is 
shown in Fig. |l(| Although the FP action is not designed to improve asymptotic 
scaling, just scaling, we can still see if asymptotic scaling obtains. Using the 
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original coupling in the action, we compute T c /A from the two loop formula from 
our data, and record the results in Table ^ and Fig. [□]. The eight parameter FP 
action not only shows asymptotic scaling within twenty per cent from Nt = 2 
to 6 but in addition the value of its A parameter is about a factor of ten larger 
than the Wilson one and therefore much closer in value to the continuum A 
parameters. 
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Figure 9: Deconfincmcnt fractions for the eight-parameter FP action for N t = 4 
lattices. 
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Table 6: Critical couplings at finite volume and extrapolated to infinite volume 
for the FP action with parameters in Table 1. 



volume 


N t = 2 


N t = 3 


N t -4 


N t = 6 


4 3 


2.89(1) 








6 3 


2.945(15) 


3.320(15) 






8 3 


2.96(1) 


3.34(1) 


3.50(1) 




10 3 


2.960(5) 


3.34(1) 


3.51(1) 


3.78(2) 


12 3 


2.975(5) 


3.343(5) 


3.520(5) 


3.775(5) 


15 3 








3.777(5) 


18 3 








3.790(5) 


infinite 


2.976(5) 


3.345(5) 


3.520(5) 


3.790(6) 


T c /A 


7.68 


7.61 


6.87 


6.13 



4.4 Torelon masses 

The string tension a is measured through the correlator of pairs of Polyakov 
loops (or more complicated objects; see below): on a lattice of transverse size L 
the correlator of two Polyakov loops averaged over transverse separations and 
separated a longitudinal distance z is 

C(z) = ^Re(P J (r ± ,z)P J t (0,0)) ~ exp(- M z), (9) 

r± 

where /1 is the so-called torelon mass. On an infinite lattice fi = oL and we will 
make the same definition on a finite lattice. 

A problem with Polyakov loop correlators is that the signal to noise ratio 
rises exponentially with the string tension, L, and z: a simple calculation shows 
that with N measurements 

Signal/Noise ~ Vn exp(-aLz). (10) 

Since as [3 falls, a rises, this means that at small f3 it is difficult to observe a 
signal unless z is very small. However, at small z the signal may be contaminated 
by excited states, and it may be impossible to go to large enough z to be certain 
that one is seeing the asymptotic behavior of the signal before it disappears in 
the noise. 

We are primarily interested in lattice spacings in the range aT c = 1/2 to 1/8 
(for the Wilson action, 5.1 < j3 < 6.0). At large lattice spacing (aT c = 1/2) 
the best signals came from the Polyakov loop itself, measured with the Parisi 
- Petronzio - Rapuano [ p"5| multihit variance reduction method. At smaller 
lattice spacing (aT c < 1/3) the signal from this operator degrades and we had 
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more success using correlators of APE-blockcd |l(| links: we iterate 

V 3 n+1 (x) = (1 - a)V?{x) + a/A^WWix + k)Vp(x+j)^ 

+ Vg{x-kjtV?{x-k)V?(x-k + j)) (11) 

(with V®{x) = Uj(x) and V J n+1 (x) projected back onto SU(3)), with a varying 
from 0.2 at aT c = 1/3 to 0.5 at aT c = 1/4 to 0.7 at aT c = 1/8 and the number 
of blocking steps rising from 10 to 15. Presumably a large part of the noise is 
local (in physical units) and hence more nonlocal (in units of a) averaging is 
needed as the lattice spacing falls. 

The multihit algorithm is not very convenient when using the FP action 
on a parallel machine for the same reason that Metropolis or overrelaxation is 
inconvenient, but as we only ran the small lattice spacing simulations on parallel 
machines that is not a problem. 

We determined the torelon mass from a single-exponential correlated fit, 
beginning at a minimum z where the \ 2 /DF is near unity, and where we see 
stability in the effective mass over a range of z. At aT c = 1/2 the fits are to 
z = 1 and z = 2 points only since there is no signal for z > 2. 

As a check of our programs we reproduced the torelon mass measurements 
of Michael and Teper on Wilson f3 = 5.9, 12 3 and (3 = 6.0, 16 3 lattices. 
At lower (3 there is little Wilson data for torelons which is a variational bound 
to which we can compare: at = 5.69 our torelon mass is greater than the 
= 5.7 result of de Forcrand et. al. (which however was done using the 
no n- variational cold wall source method.) 

An uncertainty in the value of the critical coupling propagates into the un- 
certainty of G for the FP action. By doing lower-statistics simulations at cou- 
plings close to the critical /3-values, we can roughly compute the variation of 
the torelon mass with (3. For the uncertainties in /3 c 's recorded in Table ^, this 
is an additional uncertainty in G of about 0.04-0.05 at r = 2 and r = 3/2. We 
include this uncertainty in the figures by combining it in quadrature with the 
statistical fluctuation in G; in the figures the extreme range of the vertical error 
bars shows the combined uncertainty. 

The result of our simulations is that the torelon mass measured on aspect 
ratio 2 lattices using the FP action scales within the statistical errors for 1/4 < 
a>T c < 1/2, at a value which is consistent with the value inferred from Wilson 
action results from small lattice spacing simulations. 

As a second test of scaling we compare torelon masses measured on lattices 
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Table 7: Our measurements of torelon masses from the Wilson action on small 
lattices. 



volume 


a 


fi = La 


T J 

L^ 


3 3 x 16 


5.10 


1.83(3) 


2.34(2) 


4 3 x 16 


5.10 


2.86(6) 


3.38(4) 


6 3 x 16 


5.55 


1.65(5) 


3.15(5) 


6 3 x 16 


5.69 


0.75(2) 


2.12(3) 


8 a x 16 


5.69 


1.11(2) 


2.98(3) 


9 3 x 16 


5.9 


0.45(1) 


2.03(1) 


12 3 x 16 


6.0 


0.44(1) 


2.31(4) 


Table S 


\: Our measurements of torelon masses 


from the FP action. 


volume 


/9 


fi = La 


L^ 


3 3 x 12 


2.975 


1.50(6) 


2.12(4) 


4 3 x 8 


2.975 


2.38(5) 


3.09(3) 


6 3 x 12 


3.34 


1.55(5) 


3.05(5) 


6 3 x 12 


3.52 


0.78(1) 


2.16(1) 


8 3 x 12 


3.52 


1.19(4) 


3.08(4) 


9 3 x 16 


3.79 


0.50(1) 


2.12(2) 



of constant aspect ratio 3/2. The scaling plot is shown in Fig. g. In this 
figure the diamond shows the result of our simulation with the Wilson action 
at (3 = 6.0 and the star shows an extrapolation to P c (N t = 8) = 6.06 using 
asymptotic freedom. All the simulations were performed by us. Again, the FP 
action torelon mass scales for 1/6 < aT c < 1/2, while the Wilson action exhibits 
considerable scaling violation over the same range of lattice spacings. 



4.5 The potential 

Finally, we have measured the potential from the correlation function of Polyakov 
loops. The results obtained at f} c (N t = 2) and f3 c (N t = 4) are shown in Figs. || 
and |] for the Wilson and the 8 parameter FP actions, respectively. At large 
distances the system should forget the errors related to the bad operators and 
rotational symmetry has to be recovered. For the Wilson action at Nt = 2 
strong symmetry violation is seen even at rT c = 1.5 and there is a clear scaling 
violation in the long distance part of V(r). These cut-off effects drop below the 
statistical errors for the FP action. 
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5 Conclusions 



We have described a systematic program for constructing fixed point actions 
for SU(3) gauge theory, and illustrated it using a particular scale two blocking 
transformation. An approximation to the FP action shows scaling within our 
small statistical errors beginning at aT c = 1/2, as compared to aT c < 1/8 for 
the Wilson action, at a cost of a factor of about 7 in computational speed per 
site update. 

This work allows one to consider a large number of questions. 

First, the action which is simulated here and shows scaling is an approxima- 
tion to a FP action, not to an action lying on any known renormalized trajectory 
at small (3. Can one find the renormalized trajectory of these renormalization 
group transformations? We believe that this is feasible using the available tech- 
niques |l9|. As a related question, the FP action scales for j3 values at which 
it was not a good fit to the FP action. It would be useful to understand how 
imperfect an action can be and still perform well. Additionally, it might be that 
completely different scale transformations, such as the y/3 transformation po| , 
might yield FP actions which are even more local than this one. 

In the meantime, the fixed point action can be used for many phenomeno- 
logical calculations at large lattice spacing. 

A particularly severe and unexpected bottleneck in this calculation was the 
absence of good scaling tests using the Wilson action. We are particularly 
hampered by imprecision in published values of deconfincmcnt couplings for 
N t > 8. We urge groups doing simulations with the Wilson action to choose 
parameters which are easy to relate to other length scales. We would also like to 
encourage developers of alternative improved actions to test them at very large 
values of the lattice spacing, since either scaling violations will be large and 
easy to see, or they will be small, showing that the candidate action is indeed 
an improvement. 

Finally we would like to comment on lattice perturbation theory. The con- 
struction of the FP actions we have described here does not involve lattice 
perturbation theory in the sense that we have made no expansions in g 2 or in 
powers of the lattice spacing a. The actions themselves are not explicitly de- 
signed to have good perturbative expansions, although the one we have tested 
most extensively probably does so in terms of its bare coupling (since it seems 
to scale asymptotically with a small T c /A ratio). Perturbative calculations for 
these actions are unpleasant but not impossible, given symbolic manipulation 
programs. Direct numerical calculations of renormalization constants |2l| are 



23 



in principle no more difficult than for the Wilson action, although the currents 
themselves must be "perfect." 

These actions appear to open the way towards doing QCD simulations with 
modest computer resources. 
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A Appendix 

We discuss here briefly the arguments which show that the numerical analysis 
of the classical FP equation Eq. (|J) can be done on small lattices. 

Consider a lattice of size (mL) 4 onto which a periodic L 4 coarse V con- 
figuration is copied m 4 times, where m is a large integer. The corresponding 
minimizing configuration Uo(V) lives on an (2mL) 4 lattice. The fine links U^n) 
and U ' fj,(n + 2Lv) are influenced identically by the coarse field V. Assuming that 
this symmetry is not broken spontaneously, the minimizing configuration will 
be periodic with period 2L. The relation between the coarse field and the min- 
imizing fine field on the big lattices is fixed by the relation between the L 4 and 
(2L) 4 lattice. On the other hand, the big (mL) 4 configurations are certainly 
legitimate candidates to be used to parametrize the FP action, since Eq. (||) 
should be satisfied on any type of V configurations. 

The question remains whether the specific configurations considered above 
give sufficient information to fix the coefficients in Eq. (||). In spin models it 
is easy to see that solving the FP equation on small lattices one is not able 
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to resolve certain couplings in the FP action. Working on a 2 2 lattice in the 
non-linear cr-model, for example, one can not separately determine the two- 
spin couplings p(l, 0), p(3, 0), p(5, 0), ... - only their sum enters. Using a limited 
number of loops in a non-Abelian gauge theory on a small lattice it is not 
clear whether such relations exist. We emphasize, however that even if they 
exist, no error is caused by small lattices, only information is lost. Since the 
general parametrization in terms of loops is highly redundant and our specific 
parametrization is not at all systematic, this is not a real problem. 
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Figure 10: Volume dependence ol critical couplings for the eight parameter FP 
action, and their extrapolations to infinite volume. 
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Figure 11: T c /A for the Wilson (octagons) and FP (squares) actions. 
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